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Abstract 

Quantum fields are well known to violate the weak energy condition of general 
relativity: the renormalised energy density at any given point is unbounded from 
below as a function of the quantum state. By contrast, for the scalar and elec- 
tromagnetic fields it has been shown that weighted averages of the energy density 
along timelike curves satisfy 'quantum weak energy inequalities' (QWEIs) which 
constitute lower bounds on these quantities. Previously, Dirac QWEIs have been 
obtained only for massless fields in two-dimensional spacetimes. In this paper we 
establish QWEIs for the Dirac and Major ana fields of mass m > on general 
four-dimensional globally hyperbolic spacetimes, averaging along arbitrary smooth 
timelike curves with respect to any of a large class of smooth compactly supported 
positive weights. Our proof makes essential use of the microlocal characterisation 
of the class of Hadamard states, for which the energy density may be defined by 
point-splitting. 



1 Introduction 



In general relativity, it is customary to assume that the stress-energy tensor satisfies one 
or more of the classical energy conditions; the weak energy condition, for example, being 
the assertion that the energy density measured by any observer is nonnegative. The 
primary motivation behind these energy conditions is that they ensure that gravity acts 
as an attractive force (in the sense of focussing geodesic congruences) in accordance with 
our experience of gravitation on a wide range of scales. It is therefore natural to assume 
that physically reasonable forms of classical matter obey such conditions, and to regard 



matter theories (such as the nonminimally coupled scalar field |Q, ^ |T2|, |T7l) violating 
such conditions as being of questionable physical significance on many scales. Moreover, 
the energy conditions have proved to be of great value in obtaining deep results in classical 
general relativity, such as the positive mass and singularity theorems ^ pT[ ]. 

However, it is well known that all the pointwise energy conditions are violated in 
quantum field theory. Indeed, Epstein, Glaser and Jaffe ^ proved that no Wightman 
field theory on Minkowski space can admit a (nontrivial) energy density observable whose 
expectation values are bounded from below and vanish in the Minkowski vacuum state. 
Moreover, in linear field theories (both in fiat and curved spacetimes) it is easy to con- 
struct states whose energy density at a given point may be tuned to be arbitrarily nega- 
tive P, BH]. This raises the possibility that quantum matter might be used to construct 



spacetimes with exotic properties, such as traversable wormholes [jTSf or so-called 'warp 
drive' spacetimes [Q, usually excluded by the classical energy conditions. One might 
also ask whether the conclusions of the singularity theorems remain valid for quantum 
matter. Furthermore, it is necessary to understand how classical matter contrives to obey 
the classical energy conditions, given that its fundamental constituents need not. 



One profitable line of enquiry, starting with the work of Ford IQ, has been to in- 



vestigate weighted averages of the renormalised energy density along the worldline of an 
observer, or over a small spacetime region. It turns out that the expectation values of 
these averaged observables are bounded from below independently of the state, and such 
bounds have been developed in successively greater generality over the last few years 
TB| , [TT], 1^, ^ |22|, ^. Most recently, one of us has established the existence of 



such bounds (and given an explicit, though not optimal, lower bound) for the minimally 
coupled real linear scalar field in any globally hyperbolic spacetime, in the case where 
averaging is performed with respect to proper time along any smooth timelike curve using 
an arbitrary smooth compactly supported positive weight belonging to the classQ 

W = {f eC^{R)\ fir) = g{Tf for some real- valued geC^{R)}. (1.1) 

The constraints imposed by these lower bounds have been called 'quantum inequal- 
ities' by various authors. However, as we hope to discuss elsewhere, there seem to be 
strong parallels between the phenomena discussed above and various situations arising 
in quantum mechanics. Indeed, quantum inequalities appear to be a widespread feature 
of quantum theory as a whole, stemming ultimately from the uncertainty principle. In 
this light, we adopt the more specific terminology 'quantum weak energy inequality' (or 
QWEI) in relation to lower bounds on the renormalised energy density of a quantum field. 



^See the remarks following Theorem 4.1 for a brief discussion of this class 



1 



There are also related constraints which demand that the integral of the energy density 
over any complete (i.e. inextendible) smooth timelike or lightlike ('null'-) curve be non- 
negative. These are called the 'averaged weak energy condition' (AWEC) or 'averaged null 
energy condition' (ANEC), respectively; they may be viewed as a limiting case of QWEIs 
when the weight function (/ in eqn. ( |1.2| ) below) against which the energy density is 
integrated along a timelike or lightlike curve approaches the unit function. Historically, 
it was first pointed out in a work by Tipler that suitable versions of AWEC or 
ANEC imply singularity theorems in general relativity similar to those which one obtains 
from pointwise positivity conditions on the energy density. Subsequently, the question 
of whether quantum fields obey AWEC or ANEC has been investigated in a number of 



works p6|, 0, |50|, |T^, ^ . Most of these references treat linear quantum fields; 



i2| establishes ANEC for general (axiomatic) quantum field theory in two-dimensional 



Minkowski spacetime. A study of the interrelations between averaged energy conditions 
and QWEIs is contained in [jr3|. We refer to [|T^ for further discussion and review of 
averaged energy conditions. 

QWEIs place stringent constraints on attempts to generate exotic spacetimes [|l^, |3^ 
and may open a route towards proving results analogous to the singularity theorems for 
quantum matter To date, however, most QWEI results have been obtained 



for scalar field theories, while the more physically interesting electromagnetic and Dirac 
fields have received comparatively little attention. Ford and Roman have considered the 
electromagnetic field in Minkowksi space |T6| and shown that a QWEI holds for averaging 
against the Lorentzian weight /(r) = tq/[ti{t'^ + Tq)] along timelike geodesies; this result 
was generalised to static trajectories in static spacetimes by Pfenning |]^, who has also 
removed the restriction to Lorentzian weights It is reasonable to suppose that even 
more general QWEIs may be obtained for this case. 

Both the scalar and electromagnetic fields have the property that the classical en- 
ergy density is manifestly nonnegative, a fact which underpins all the results on these 
fields. The Dirac field is technically very different in that the 'classical' energy density 
is unbounded both from above and below; in second quantization, renormalisation serves 
the dual purpose of restoring finiteness and imposing positivity of the Hamiltonian. This 
problem appears to have restricted progress on the Dirac field to date. The main contribu- 
tion has been that of Vollick, who established a QWEI for Dirac fields in two-dimensional 
spacetimes ES] by converting the problem to one involving a scalar field and then adapting 



arguments due to Flanagan [|TT| . There seems little prospect of generalising this argument 



beyond the two-dimensional setting. In four dimensions, Vollick has also given explicit 
examples of states with locally negative energy densities |0 and demonstrated that the 
resulting energy densities nonetheless obey QWEIs modelled on those for the scalar field. 

In this paper we establish a general QWEI for massive or massless Dirac fields on 
four-dimensional0 globally hyperbolic spacetimes. To be more specific, let 7 be a smooth 
timelike curve, parametrized by its proper time r, in a globally hyperbolic spacetime 
{M,g). Let ujq be a given (but arbitrary) Hadamard stateQ of the Dirac field on (M, gr). 
The state uq is used as a 'reference state' to define the expected normal ordered energy 

^The restriction to four dimensions is purely for convenience: our methods would apply in more general 
dimensions. 

■^See Sec. || for a brief review of the concepts used here. 
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density (: Tqo ■)uj for any other Hadamard state uo. Our main result, Theorem asserts 
that 

inf j dr (: Too :)^(7(r))/(r) > -oo, (1.2) 
where the infimum is taken over the class of Hadamard states and / belongs to the class 



W. In principle our arguments yield an explicit lower bound for the left-hand side of ( p.. 21) . 
This expression is unfortunately not particularly enlightening and is not expected to be 
sharp. Let us note that ( |1.2|) remains true if the normal ordered energy density is replaced 
by the renormalised energy density, as these two quantities differ by a smooth function. 

The plan of the paper is as follows. We begin, for completeness and to fix notation, by 
reviewing the theory of the quantized Dirac field in Sec. 0. Particular attention is given to 
the class of Hadamard states, which may be characterised by a microfocal spectrum con- 
dition on the wave-front set of the two-point function. This formulation of the Hadamard 
condition is technically convenient and allows us to bring the tools of microfocal analysis 
to bear. Sec. |] explains how the normal ordered energy density may be constructed by 
point-splitting. 

The proof of our QWEI begins in Sec. using the following strategy. The averaged 
normal ordered energy density is first expressed as an integral over R^; decomposing this 
integral according to the quadrants of M^, each piece is then split further into four using 
a decomposition induced by the reference state uq. All but two of the resulting sixteen 
contributions can be bounded (both above and below) using estimates obtained in Sec. |^. 
The remaining terms are then expressed in the form TZ = limA_^_|_oo Tr JaW where Ja and 
W are self-adjoint and Ja is independent of lj. The parameter A G M"*" defines a cut-off, 
used to avoid domain problems. We prove that W is positive and trace-class with bounded 
trace as u varies. To conclude that TZ is bounded below, it then suffices to establish that 
the operators Ja are bounded below uniformly in A. This is accomplished in Sec. 
completing the proof of our QWEI. In Sec. ^ we briefly describe how our arguments can 
also be applied to the Majorana field. 

Conventions: The metric signature is (+,—,—,—). Lower (resp. upper) case Latin 
characters from the beginning of the alphabet will label tetrad (resp. spinor) indices. 
Tetrad indices run from to 3, and we will use j, k to label the spatial components 1, 2, 3. 
The summation convention will be used throughout the paper except where otherwise 
indicated. Units with c = h = 1 are adopted. 

The Fourier transform of an integrable function / on M" will be defined using the non- 
standard convention f{k) = J d"'x e^''^ f{x), with inverse h{x) = (27r)~" J d^ke~'^^^h{k). 
The Fourier transform of a distribution u with compact support is u{k) = u{ek) with 
ek{x) = e^^^. 

Given a Lorentzian manifold {M,g), S)'{M) will denote the space of distributions on 
M as defined in §6.3 of [^. Thus if m G ^'(M) there is, for each chart (f/, k) in M, a 
distribution G ^'{n{U)) such thatQ 

n{f) = u^ii^fMf) o K-') yfeC^iU) (1.3) 



*The factor of y— jsl used to identify test functions with test densities, on which u strictly speaking 
acts. 
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and so that = {k' o k,~^)*u'^ in k,{U fl U') for any other chart {U', k,'). We will also 
(and usually) write u o k,~^ for m^. Spinor and cospinor distributions will be defined in 
an analogous fashion, with the convention that a spinor distribution acts on test cospinor 
fields, and vice versa. 

2 The Quantized Dirac Field 
2.1 Geometrical preliminaries 

In order to make the present paper sufficiently self-contained, we need to summarize a 
few basic facts about the geometry of spinor fields in curved spacetimes. We will follow 
Dimock's work to large extent. 

We will consider Dirac fields in a four- dimensional globally hyperbolic spacetime 
{M,g). To begin with, we recall that a globally hyperbolic spacetime is a Lorentzian 
spacetime admitting a Cauchy surface, the latter being a smooth hypersurface in M 
which is intersected exactly once by each inextendible causal curve in (M, g) . We will 
also suppose that (M, g) is orientable and time-orientable, and that such orientations 
have been chosen. Then (M, g) possesses a spin-structure, that is, there is a principal 
fibre bundle S{M,g) having SL(2, C) as structure group, acting from the right, together 
with a 2-1 fibre-bundle homomorphism : S{M,g) —>■ F{M,g) which projects S{M,g) 
onto the frame bundle F{M,g). That is to say, preserves base-points and obeys 

^Oi?, = i?A(s)0V', (2.1) 

where R denotes the right action of the structure groups on the principal fibre bundles 
involved and SL(2,C) 9 s i— >■ A(s) G =SfJ is the covering projection onto the proper 
orthochronous Lorentz group. We recall that F{M, g) is the bundle of oriented and time- 
oriented tetrads (cq, ei, 62, 63), so that g{ea, e^) = rjab with rjah = diag(l, —1, —1, —1) where 
Co is timelike and future-pointing and the tetrad is given the orientation of M. 

Moreover, a collection of 4 x 4-matrices 70, . . . , 73 is called a set of Dirac matrices if 

7a7f. + 7fe7a = 2?7afe " 1 . (2.2) 

A theorem due to Pauli states that, if 70, • • • ,73 and 7o, • • • ,73 are two sets of Dirac 
matrices, then there is an invertible matrix M so that 7a = M7^M~^. Any set of Dirac 

A( . ) t 

matrices 70, • • • ,73 is connected to the covering SL(2, C) ^ in the following way. 
Let Spin(l, 3) consist of all unimodular 4 x 4-matrices S so that 

SjaS-' = (2.3) 

holds for some real numbers A^^ = ^^a{S)- It follows from the defining properties of 
Dirac matrices that K^aiS) is contained in the Lorentz group. The restriction of the map 
S ^ A.^a{S) in ( p.3| ) to SpinQ(l, 3), the unit connected component of Spin(l, 3), is a group 
homomorphism with range .ifj, and thus Spino(l,3) is isomorphic to SL(2,C). 

Sometimes it is useful to distinguish sets of Dirac matrices with certain properties. 
One says that a set 70, . . . , 73 of Dirac matrices belongs to a standard representation if 

7o = 7o and 7^ = -7fc . (2.4) 
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Here, 7* is the hermitean adjoint of 7a. A set of Dirac matrices which belongs to a standard 
representation and has the additional property that the complex conjugate matrices fulfill 

7„ = -7a, (2.5) 

is said to belong to a Majorana representation. 

We will now suppose that we are given a globally hyperbolic spacetime (M, g) together 
with a spin-structure {S{M,g),ip) and a set of Dirac matrices 70, • • • ,73 which we will 
assume, for the sake of notational simplicity, to belong to a standard representation. (We 
note, however, that everything which follows could also be carried out in a similar way 
without that assumption.) As was pointed out above, via the isomorphism SpinQ(l, 3) ~ 
SL(2, C), carries a representation of the universal covering group of J/f^. which is given 
by the action of the matrices S in SpinQ(l, 3) on vectors in C"^. Via Spino(l, 3) ~ SL(2, C), 
we can also regard S{M, g) as a SpinQ(l, 3)-principal bundle and form the associated vector 
bundle 

DM = S{M,g) ><spino(i,3) (2.6) 
That is, the fibre of DM at p e M consists of the orbits 

[Sp,x] = {{Rs%,Sx) : S e Spino(l,3)} (2.7) 

for Sp G S{M.g)p and x G C^. There is a fibrewise left action of SpinQ(l,3) on DM by 
Ls[s:p^x\ = [Sp, Sx] . Elements in DM are called spinors, and elements in the dual bundle 
D*M are called cospinors. Moreover, if is a (local) section in S{M,g), then it induces 
on one hand a tetrad field (cq, ... ,63) = ip o E, i.e. a (local) smooth section in F{M,g), 
via the spin-structure, and on the other hand it induces a set {Ea)a=i of (local) smooth 
sections in DM, defined by 

Ea = [E, 6a] (2.8) 

where bi, . . . ,64 is the standard basis in C^. There are corresponding dual tetrad fields 
(e°,... ,6^) defined by e^(e„) = and {E^)%=i defined by E^{Ea) = S^. The are 
smooth sections in T*M, and the E^ are smoth sections in D*M, the dual bundle to 
DM. We shall denote by C°°{DM) and C°°{D*M) the sets of smooth sections in DM 
and D*M, respectively. The notation for smooth sections in TM and T*M will be similar. 

With respect to the given set of Dirac matrices, one can define a section 7 in C°°(T*M)® 
C°°{DM) ® C°°{D*M), i.e. a mixed spinor-tensor field, by setting its components •jh^B 
in the induced frame ® Ea'^ E^ to be equal to the matrix elements {'yb)^B of 7^. This 
definition is independent of the induced frames, i.e. independent of the chosen (local) 
section E in S{M, g). (Once the set of Dirac matrices 70, ■ ■ ■ ,73 is given, 7 encodes the 
spin-structure at the level of DM.) 

Moreover, there is an anti-linear isomorphism DM — > D*M induced by forming the 
Dirac adjoint: U u — u^Ea is a spinor, one can assign to it a cospinor — u^E^ with 
components 

it+ = M^7oAB (2.9) 
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where Jqab are the matrix elements of 70. This assignment possesses an inverse (denoted 
by the same symbol), where a cospinor v — vbE^ is mapped to a spinor v+ having 
components f = 'Jq^^vb- Again, jo^^ are the matrix elements of 70. The operation 
of taking the Dirac adjoint gives rise to anti- linear isomorphisms between C°°{DM) and 
C'^{D*M) in the obvious manner. 

The metric- induced covariant derivative V on C°°{TM) induces a covariant derivative, 
also denoted by V, on C°°{DM). If (eo, ... ,63) and {Ea)a=i are induced by a section E 
in S{M,g) and / = J'^Ea is a local section in DM, then 

V/ = V6/^(e^ ® Ea) e C~(r*M) ® C°^{DM) (2.10) 

has the frame components 

VJ'" = 9,r + a/B/^, where dj"" = df^'ie,) , a/^ = -^r«,7/c7% • (2.11) 

Here, df^ is the differential of the function f^, and we read the components of the Dirac 
matrices on the right hand side while F^^ are Christoffel's connection coefficients, defined 

by 

VA; = (a^r + T'^dky ® ea (2.12) 

for k = k^eb G C°^{TM). 

The covariant derivative V can be extended to cospinor fields and mixed spinor-tcnsor 
fields by requiring the Leibniz rule and commutativity with contractions. Thus, ii h = 
HbE^ is a cospinor field, then Vh — Vbhse^ ® E^ has the components 

VbhB = dbhB - hcab^B ■ (2.13) 

It follows that V7 = 0. 

2.2 The Dirac Equation 

The Dirac-operator ^ is a first order differential operator taking spinor fields to spinor 
fields, or cospinor fields to cospinor fields; it is defined as the action of the covariant 
derivative followed by contraction with the spinor-tensor 7. More precisely, ii f — J^Ea £ 
C°°{DM) and h = HbE^ e C°°{D*M), then 

y/ = (m^EA^V^'la^B^bf^EA, (2.14) 
= {fh)BE^ ^ri'^Vbhc^a^BE^ . (2.15) 

An important property of the Dirac operator is that it commutes with taking the Dirac 
adjoint, i.e. 

(y/)+ = y/+ and {fh)+^fh+ (2.16) 
for / e C°°{DM) and h e C°°{D*M) . 
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The Dirac equation is the following first order partial differential equation for spinor 
fields u G C°°{DM) or for cospinor fields v E C°°{D*M): 



-i^ + m)u 



0, 
0, 



where m > is a constant. Then 

P 



m)(-iy+ m) 



(2.17) 
(2.18) 



(2.19) 



is the Lichnerowicz wave operator on spinors or cospinors. It is a second order wave 
operator which has metric principal part, and owing to global hyperbolicity of (M, g) this 
implies that the Cauchy problem for the corresponding wave equations is well-posed and 
that P possesses uniquely determined advanced and retarded fundamental solutions. As 
shown in 0, this implies that the Dirac operators —i'^+m on spinor fields and + m 
on cospinor fields possess uniquely determined pairs of advanced(— ) and retarded(+) 
fundamental solutions and respectively: This means that, for the spinor case. 



sp 



C^{DM) 



are continuous linear maps so that 

(— m)Sfu 



u 



C^{DM) 



(2.20) 



(2.21) 



holds for all u G C^{DM) and, moreover, supp S^u is contained in the causal fu- 
ture (+) /causal past(— ) of suppw. (We note that our convention concerning advanced/re- 
tarded fundamental solutions is opposite to that in ||3^.) The cospinor case is analogous. 
Then one defines the retarded-minus-advanced fundamental solutions 



5'sp = — and S, 



sp 



■'cosp 



q+ _ C- 

"-"^cosp "-^cosp 



(2.22) 



In order to quantize the Dirac field, it is very convenient to 'double' the system by 
taking pairs of spinor fields and cospinor fields together, as was done in the references 
We give here an equivalent version which makes contact with the notation 



28 



23 



used in [g. To this end, let us denote C^{DM) by %p and C^{D*M) by ^cosp, and 

© %p. On ^double we introduce the sesquilinear 



define the doubled space ^double 
form 



^cosp 









' h2 ' 


( 
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{ftJ2)-{hi,ht) 



(2.23) 



for hi,h2 G 
pairing 



^cosp and /i, /2 G ^sp, where for v G ^cosp and u G ^sp we employ the dual 



{v,u) 



(2.24) 



M 



with dfig denoting the canonical 4-volume form induced by the metric g on M. This 
dual pairing also embeds ^cosp in the (topological) dual space of ^sp, and vice versa. 
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embeds in ^cosp- The sesquilinear form ( . , . ) is non-degenerate, but not positive. A 
useful relation is 



(2.25) 



for h e ^cosp and f e %p. 

Let us define the conjugate-linear isomorphism T : ^double ^double, playing the role 
of a charge-conjugation, by 



(2.26) 





' h ' 






( 


_ f . 


)- 





Then one finds that F is a skew-conjugation with respect to ( . , . ' 

(m, TFs) = -(F2, Fi) ^F,,F2e ^double 

Now we introduce the following 'doubled' operators on ^double: 



namely it holds that 
(2.27) 



-y+im 
y+ 



—'^—im 
im 



^cosp 

iS. 



sp 



Then it holds that 



D<iD^ — D^D^ — Pdouble 



P 

cosp " 

R 



sp 



(2.28) 
(2.29) 

(2.30) 



where P., denotes the Lichnerowicz wave operator on spinors and cospinors, respectively; 
moreover, one finds that F commutes with Pdoubie and 



(2.311 



-(Pi, D^F2) , {D^, Pi, P2) = -(Pi, P^Ps) \fFi,F2e ^double . (2.32) 



One may also check that (defined in obvious analogy to S^) are the retarded(+)/ad- 
vanced(— ) fundamental solutions for the operator Pt>; consequently 

Furthermore, from (|2.31|) one can see that 

FS< 

This entails that F is a complex conjugation for the sesquilinear form 

(Pi, F2)s = ('S'<]Pi, P2) , Pi, P2 G ^double , 

SO that 

(FPi,FP2)5 = (P2,Pl)5 ^F,,F2e 

■^double ■ 



(2.33) 



(2.34) 



(2.35) 



(2.36) 
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On the other hand one can see that (cf. ( p.25| )) 



( ^1 , ) =-^{ft,SsJ2)+^{Scos:>h,ht) (2.37) 

and this imphes by Prop. 2.2 in ||^ that ( . , . )5 is positive-semidefinite, {F,F)s > 0. 

Now we introduce the quotient space ^double /ker and denote by its comple- 
tion with respect to ( . , . )5. The conjugation T induces by (|2.34|) a conjugation on 
which will again be denoted by F. Hence, we have derived from the doubled Dirac equa- 
tion a complex Hilbert-space Jif (with scalar product ( . , .)s) together with a complex 
conjugation F. The system can be quantized, following Araki [|l|, by assigning to these 
data the algebra of canonical anti-commutation relations CAR.{Jif, F). This is the unique 
C*-algebra with unit 1 which is generated by a family {-B(v) : v G Jif} subject to the 
relations: 



(i) V I— s> -B(v) is C-linear, 

(ii) i?(Fv) = Biv)* , 

(iii) B{y)*B{w) + B{w)B{^r)* = (v, w)^ ■ U . 



/e^sp. (2.39) 



Now let q : ^double ^doubie/ker denote the quotient map, then we define the quan- 
tized Dirac field as the linear map which assigns to each h G ^cosp the element 

^(/^) = i?(?[o ]) 

in CAR{Jf, F). The adjoint spinor field will be defined by 



(2.38) 





' " 




/ 



As a consequence of (iii), the field and its adjoint satisfy the anti-commutation relations 

^{h)^^if) + ^^imih) = -t{h, 5sp/) . (2.40) 
We also note that, owing to ( p.33| ), B{q{D^F)) = 0, and this entails 

*((iy+m)/i) = and ^+((-iy + m)/) = (2.41) 
for all h G ^cosp and / G 

Remarks, (i) acts linearly on cospinors and fulfills, in the sense of distributions, the 
equation (— m)\E' = 0, therefore the map h i— >■ "^{h) is regarded as a spinor field. 
Similarly, \E'"'" acts linearly on spinor fields and fulfills in distributional sense (i'^+m)'^'^ = 
0; hence it is viewed as a cospinor field. 

(ii) ^ and \E^+ are C*-valued distributions since e.g. 2||\&+(/)|p = -«(/"*", 5'sp/) and 
%p ® ^sp 3/1^/2'-^ ~^{fii 'S'sp/2) is continuous (with respect to the usual test-function 
topology). 
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(iii) We briefly comment on the case where 70, • • • ,73 belong to a Majorana representation 
(as considered in [^]) and one wishes to quantize the Majorana field. In that situation, 
^cosp (and similarly, ^gp) carries an 'intrinsic' charge conjugation F given hj v = vbE^ h-> 
vbE^ in any frame. Then F is a skew-conjugation for the sesquilinear form 



(K h') 



dfigip) h'p{h+) 



(2.42) 



M 



on ^cosp- Upon defining = ^—im, = '^+im and = iScosp, one obtains similar 
relations as before. Then arises as completion of ^cosp/ker and {h, h')s = {S^^h, h'). 
The field operators then simplify to 



^>{h) = B{q{h))- vI/+(/) = vI/(/H 



v^(F/+) 



(2.43) 



-'cosp ) 



/ G ^sp- In other words, the Majorana field may be quantized without 



for h G 

doubling the classical system. Instead of starting with ^cosp one can likewise consider 
S>sp'-i this has been done in ^ . 



2.3 Hadamard states 

We recall that a state on a C*-algebra C is a linear functional uo : C ^ £, fulfilling uj{t) = 1 
and u!{A*A) > for all A & C. For the purpose of the present work, it is sufficient to 
focus on the two-point functions 002 of states u on CAR(t^, F). The two-point function 
t02 of to is an element in (i^doubic ® -^double)' given by 



tL!{B{q{Fi))B{q{F2))) , -Fi,-F2G ^double • 



(2.44) 



It was shown in |l|] that there is for each state uj on CAIi{Jif, F) a linear operator Q on 
with the properties: 

(I) < g* = g < 1 , 

(II) g + FgF = 1 , 

(III) uj2{Fi ® F2) = (Fg(Fi), Qq{F2))s for all Fi, F2 G ^double- 

Conversely, each linear operator Q with the properties (I) and (II) determines by (III) 
the two-point function 002 of some state u on CAR(^, F), a so-called quasifree state, 
determined by the two-point function 002 (see [|l| for discussion). Such a quasifree state u 
is pure (and then often called a Fock-state) if and only if g is a projection, i.e. = Q. 
The operator Q with the properties (I), (II) and (III) above will be called the operator 
labelling the quasifree state u. 

The stress-energy tensor is defined using the two-point functions of a particular class 
of states, the Hadamard states. One says that a state u on CAR(J^, F) is a Hadamard 
state if we may write 



uj2iFi ® F2) = w{D^Fi ® F2) 



(2.45) 
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for some distribution w G (^double ® -^double)' of Hadamard form for the doubled wave- 
operator Pdoubio on ^double! the definition of a Hadamard form for such a wave-operator 
has been given in |^ (cf. also |2^, ^). This definition entails that the difference 



between the two-point functions of two Hadamard states is smooth. For the purposes of 
this work, we will also need the characterization of Hadamard states in terms of properties 



of the wave-front set WF(ci;2) that appears in p9|, ^ |3^, following a line of argument 



given in ||35| for the scalar case. The relevant statement, proven in the references just 
stated, is: A state uo on CAR(^, F) is a Hadamard state if and only if the wave-front set 
of its two-point function uj2 satisfies the relation 

WF (072) = , -n e T*{M X M) I ~ i e AT+j , (2.46) 

where T*{M x M) is the cotangent bundle over M x M without the zero-section, {p,C,) ~ 
{p', C,') means that there is a lightlike geodesic connecting the points p and p' in M and 
to which ^ and ^' are co-tangent, and is the set of all future-directed null covectors 
at p. (We remark that in |^ the opposite sign convention for the Fourier transform was 
chosen, leading to the opposite constraint ^ G (i.e., past-directed null covectors) in 
that reference compared to ( p.46| ).) 

At this point we very briefly recall the definition of the wave-front set of a distri- 



bution H. For a distribution t G ^'(M"), a point {x,k) e W x (]R"\{0}) is called a 
regular directed point of t if there exists x ^ ^(M") with x{^) 7^ S'-iid a conic open 
neighbourhood C of /c in R"\{0} such that 

sup (1 + |A;'|)^|£t(A;')| < oo VA^gN. (2.47) 

(If this holds we will say that xt is of rapid decay in C.) The complement in M" x (R"\{0}) 
of the set of all regular directed points of t is called the wave-front set WF (t) of t. Given 
a scalar distribution r on a manifold M, one says that a non-zero covector {p,C,) G T*{M) 
is in WF (r) if there is a coordinate chart {U, k) around p G M so that {n{p) ^) G 
WF (roK^i) where (as discussed at the end of Sec. |I|) tok^^ G S>\h(U)) is a distribution 
on the chart range of k. This definition of WF (r) is independent of the choice of the 
chart K. We refer the reader to for further discussion of the properties of wave-front 
set of distributions on manifolds. For the case that r is a distribution on test-sections of 
a vector bundle, e.g. defined on ^double, t can be viewed, via (local) trivializations of the 
bundle, as collection {t'^b)a,b of scalar distributions, and then WF (r) is defined as the 
union of WF (t'^b) over all components A^B. It is not difficult to see that this definition 
is independent of the chosen (local) trivialization. 

Now let UJ2 be the two-point function of a Hadamard state uo on GPJii{^, F), and let 
Q be the corresponding operator on M' with the properties (I),(II) and (HI) above. We 
will use the notation 

= TQT (2.48) 

for the 'charge conjugate' of Q, and we will adopt this notation also for other operators 
on In studying the stress-energy tensor, we will be particularly interested in the 
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following distributions on ® -^cosp associated with ui (respectively, with Q) , which we 
will also refer to as two-point functions: 



UQif (^h)= a;(xI/+(/)vI/(/i)) , uUf ^ h) = u;{^{h)^+{f)) 



(2.49) 



where / G and h G ^cosp- Note that Uq 



UJ, 



qf. As a consequence of the constraint 
on the wave-front set (|2.46| ) for Hadamard states, which is also called microlocal spectrum 
condition, one finds that the following microlocal spectrum condition holds for uq and 



a;, 



Q' 



WF (Jq) = {(p,e;p', -O G t*(M X M) I (p,0 ~ ip',n e e AA/} . 



(2.50) 



Here, and below, we use jl and b to denote either the presence or absence of a F in the 
following context-dependent way: 



Q UQ Af+ M+ . 

To avoid confusion we will sometimes use ■ as a placeholder to indicate the absence of a 
F. Thus Q' = Q, for example. We note that the microlocal spectrum condition in the 
form (|2.50|) has been proved directly for quasifree Hadamard states in p9|, ^ . 
Let us note that ( |2.49| ) may be written 



LUQif (^h)={Tq 




/ 



,Qq 



h 





(2.51) 



as a slight abuse of notation we will use this relation to define ujq for general bounded 
operators Q on and refer to ujq as the two-point function labelled by Q. (Of course ojq 
is not in general the two-point function of a state.) We will also denote by Had (^, F) 
the class of operators which obey properties (I) and (II) and such that uOq obeys ( ^^.50] ). 
Thus Had (J^, F) parametrizes the quasifree Hadamard states on CAR(^, F). 

Our proof of Thm. 4.1 below relies on the existence of pure, quasifree Hadamard 
states for the Dirac field. It seems that this has never been established in the literature in 
full detail, therefore we sketch here how such states may be constructed by adapting an 



argument employed by Fulling, Narcowich and Wald |jT9[ for the case of the free scalar field 
to the Dirac field. The first step is to show that there exists a pure, quasifree Hadamard 
state for the Dirac field for ultrastatic {M,g). In fact, if {M,g) is ultrastatic, then it 
may be endowed with a suitable spin structure so that the ultrastatic time shifts give 
rise to a continuous unitary group Ut (t G M) on J^f leaving the scalar product ( . , .)s 
invariant and fulfilling TUt = UtT. Moreover, if the mass parameter m appearing in the 
Dirac equation is strictly positive, then the spectrum of the self-adjoint generator of Ut 
is bounded away from zero. Theorem 2 in |jl|, or the results in |4^, thus show that 
there is a pure quasifree state ujq on the CAR algebra of the Dirac field on ultrastatic 
{M,g) which is a ground state for the C*-dynamics induced by Uf, the projection Pq 
labelling c^o is the projection onto the positive spectral subspace of the unitary group Ut- 
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Since Uq is a ground state, it fulfills the microlocal spectrum condition [0 and hence is 



a Hadamard state In a second step, one uses a technique developed in |jT9| which 
allows one to view a neighbourhood of a Cauchy surface of any given globally hyperbolic 
spacetime as being isometrically embedded in a globally hyperbolic spacetime that has 
an ultrastatic part (with suitable spin structure as above) in its past. By the uniqueness 
of the Cauchy problem and the "propagation of Hadamard form" under the dynamics of 
the Dirac equation [^, any pure, quasifree Hadamard state prescribed on the ultrastatic 
part of the spacetime (e.g. uq) induces a pure, quasifree Hadamard state everywhere on 
the spacetime, in particular on the embedded neighbourhood of the Cauchy surface of 
the initially given globally hyperbolic spacetime. Using the same argument once more, 
a pure, quasifree Hadamard state for the Dirac field is thereby induced on any given 
globally hyperbolic spacetime. The mass parameter m may be allowed to be variable over 
spacetime in this process without affecting the Hadamard form, so that one obtains a 
pure, quasifree Hadamard state of the Dirac field on any globally hyperbolic spacetime 
for any m > 0. The argument just sketched implicitly also shows that there exists an 
abundance of quasifree Hadamard states. 



3 A point-split energy density 

For the remainder of this paper, we will assume that (M, g) is globally hyperbolic, ori- 
entable and time orientable, with spin structure {S{M, g),ip) and that the Dirac matrices 
7a belong to a standard representation. 

Let 7 : M — ^ M be a smooth timelike curve in {M,g), parametrized by its proper 
time, along which we wish to establish a QWEI. The starting point is the construction of 
a normal ordered energy density on 7, which is accomplished as follows. We first claim 
that there exists a tubular neighbourhood of 7 and a local section E of S{M, g) over 

such that the induced tetrad field (cq, . . . ,63) = ip o E satisfies CqI^ = u, where m = 7 
is the velocity of 7. To see this, choose any locally finite open cover {Uj | j G Z} of 7 by 
charts Uj such that 

(i) . Uj n f/fc = unless \j - k\ < 1; 

(ii) . Uj n f/j+1 is contractable for each j G Z; 

(iii) . Uj n t/fc n f// = if j, k, I are distinct. 

The existence of such a cover follows from global hyperbolicity of (M, g) since 7 is timelike. 
Now extend m to a smooth timelike unit vector field u on some tubular neighbourhood 
of 7, so that C [JjUj. Choose any tetrad (cq, . . . ,63) on C^. Then we may 

obtain a tetrad (cq, ... ,63) with cq = m by applying a unique boost in =SfJ at each point 
(whose parameters are given by the components of u with respect to e^, and therefore 
vary smoothly). This tetrad lifts smoothly to S{M, g) in each Ujf\C^ and may be patched 
together along to obtain the required section E by virtue of properties (i), (ii) and (iii).^ 

^Of course, there are exactly two such sections. 
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Next, we choose smooth spinor fields (A = 1, . . . , 4) in C^, such that 

5^V®t^^=7o- (3.1) 

This is easily satisfied by taking va = E^, where Ea is the spin frame induced by E; 
however, it will be convenient to make a slightly different choice when considering the 
Majorana field. The changes relevant for Majorana fields will be described in Sec. |^. 
With respect to the frame Ca the Dirac stress-energy tensor is 

Tab = {^^l(aVb)i^ - (V(„V'"')76)^) , (3.2) 

which is manifestly symmetric, and is conserved provided ip obeys the Dirac equation (|2.17|) . 
In particular, Too(7('r)) is the energy density measured by an observer with worldline 7 
at proper time r. 

We may use ( p.l| ) to define a bi-scalar point-split energy density 

r{x,y) = 5^^'- {{^+VA){x){v+eo ■ V^){y) - (h ■ V^+]vA){x){v+^){y)) (3.3) 

with the property that T{x,x) = Too(x). Integrating by parts, T becomes a scalar bi- 
distribution T G (^(M) ® ^(M))' 

r(/ ^g) = 5^^'- (^+(V ■ [eoVAfmv+g) - ^P^{vAmV ■ [eov+g])) ■ (3-4) 

Here (and below) the notation V ■ [cqv] denotes minus the distributional dual of Cq ■ V, 
applied to the test function or (co)spinor v. Thus, 

V ■ [cqv] = fV ■ Co + Co ■ Vf (3.5) 

where, with respect to local coordinates (x^), V ■ Cq = V^Cq and Cq ■ Vw = CqV^w. Upon 
quantization, we obtain the algebra-valued bi-distribution T given by 

T{f 0g)= 5^^'- (vI/+(V ■ [eoVAf])^{vlg) - ^^{vA^iV ■ h^^?])) • (3.6) 

Given a state u we will now use the same symbol to denote its two-point function 
^{f®9) = '^(^^(/)^(fl')) S'lid also set vab = va®v~^. Thus vab^^ will denote the matrix 
of scalar bi-distributions 

VABU^if ^g)= u;{^^{vAm{v^g)) . (3.7) 



The formulae V-(eofA/) = f aV- (cq/) +o-o'^a/^^c and V- (eof^S') = v'I^V-{eog)+ao'^BgVc 
now allow us to write the expectation value {T)^ of T in state u as 

(t)^ = /:^Vb^, (3.8) 

where 

= 1 (1 ® ,eo ■ V - zeo ■ V ® 1) + ^6^^ , (3.9) 
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and 



qAb ^ ^[sCB^^A^ ® 1 - 1 ® 5^^ao^c] • (3.10) 

If a reference Hadamard state ujq is now specified, we may define the normal ordered 
point-split energy density (with respect to Uq) by 

{■.T% = {T)^-{TU. (3.11) 

This may also be written 

(: T :)^ = C^^'vAB : ^ : , (3.12) 

where : u := uj — ujq is the normal ordered two-point function. Since : : is smooth 
for Hadamard uj, (: T :)^ is also smooth. Accordingly, the 'coincidence limit' (i.e., the 
restriction of (: T :)^ to the diagonal) is well defined and yields the normal ordered energy 
density near 7. We denote the energy density along 7 by 

(:p:).(r) = (:T:).(7(r),7(r)). (3.13) 

Note also that (: T :)(^(/ (8> g) is symmetric in / and g by virtue of the CAR's. 

It will be convenient to regard (: p as the diagonal of the pull-back 72 (: T :)^ where 
72(r, r') = (7(r), 7(r')). In turn, 72 (: T :)^ may be written as the action of a differential 
operator on the pulled-back normal ordered two-point function: 

72*(: T:)^ = ^[{t0D-D® U)<5^^ + 72*8^^] ^^Vab : 00 : , (3.14) 

where D = id/dr (strictly speaking, D should be regarded as the distributional dual of 
—id/dr). 



4 Main argument 

We now come to the proof of the QWEI for Dirac fields. In the following, (M, g) is 
assumed to satisfy the hypotheses stated at the beginning of the previous section. 

Theorem 4.1 Let 7 : M M he a smooth timelike curve in {M,g) parametrized by its 
proper time. Let uq be a Hadamard state of the Dirac field on (M, g). Define the normal 
ordered energy density (: p by ^3.1!^ with respect to the reference state uq. Then for 
any weight f belonging to W {defined in Eq. 

inf j dT{: p:)^{t) fir) > -00, (4.1) 

where the infimum is taken over all Hadamard states uj. That is, there exists a quantum 
weak energy inequality for the Dirac field. 
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Remarks, (i) If the reference state uq is changed, (: p is modified by a smooth 
function which is independent of u. Thus we may assume without loss of generahty that 
Uq is pure and quasifree. Exactly the same argument entails that ( |4.1j ) holds if we replace 
the normal ordered energy density by the renormalised energy density. 

(ii) Perhaps surprisingly, the class W (of squares of real- valued C^(M) functions) does 
not coincide with the class of nonnegative smooth compactly supported functions. In 
fact, Glaeser ||20t| p| has constructed an example of a C°° nonnegative function /, vanishing 
only at the origin, so that ^ a//(x) diverg 0. The delicacy of this point resides 

in the behaviour of / at zeros of infinite order. It is not clear whether the restriction to 
weights in W is purely a technical limitation of our proof, or whether QWEIs should be 
understood as quadratic form results (cf. ||T0|1 ). 

Proof: It is sufficient to prove ( [4.1D for arbitrary / G C^{I) fl W where / C M is an 
arbitrary open interval with compact closure. To this end, choose r] G C^{M) such that 
r] equals unity on a neighbourhood of 7(/). It is easy to see that (: p :)^ is unaltered on / 
if we replace vab : uj : in (|3.14| ) by the compactly supported distribution uab '■ uj :, where 



uab = VVa ® V^B ■ (4-2) 

Applying the formula 

J dr F{t, tMt) = J ^H-X, X'MX - A') , (4.3) 
which is valid for F G C^(R^) and G C^(M), one may show that 

1 = J dT{:p:)M)f{r) = J dXdX' J^"" {X, X')W^/^{X, X') (4.4) 

where 

A') = [i;uAB : 00 -.f (-A, A') , (4.5) 



and we have also written 

1 

where 



■^^""(A, A') = ^ (A + A')/(A - X')5^^ + r^/]^(A - A') , (4.6) 



r^(r) = 72*0^^(T,r) = z 5^^ao^c|,M - ^^Vo^cl.M (4.7) 



is clearly hermitian (6'^^(r) = 6'^'^(r)). It follows that J^^ is a hermitian matrix kernel, 
i.e., J^^(A,A') = J^^(A',A). 

Note that J"^^ is state- independent, while contains all the dependence on the 

state of interest u and the reference state uq. We also note that J^'^(A, A') decays rapidly 
away from the diagonal in R^. 



^We are grateful to S.P. Eveson and P.J. Bushell for help in locating this reference. 
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Assuming without loss that ooq is pure and quasifree, it must be labelled by some 
projection P on J^. Since the stress-energy tensor is defined in terms of the two-point 
function, it is enough to establish ( |4.1| ) when the infimum is taken over quasifree Hadamard 
states, whose two-point functions are of the form ujq for Q G Had [J^, F) as discussed in 
Sec. The normal ordered two-point function : uq := ujq — up is labelled hj Q — P, 
i.e., : LJQ := uJq-p in the spirit of the remarks following (|2.51|) . Now 



Q-P = -PQ^P + P^QP + PQP^ + P^QP^ 
and this induces a decomposition of : ujq : into four pieces 

.p. p.. . .r r -V 

: u!q := —u! + uj +00 +00 , 
where J^^ is the two-point function labelled by P'^Q^P^] that is. 




/ 



P^Q^P^q 



h 




(4.9) 



(4.10) 



for / G %p and h G ^cosp- In Sec. ^ we will show that each cu"^^ may be pulled back to 
by 72, allowing us to write 



W 



AB 



-Wab + R + + W, 



AB 



AB 



AB 5 



where 



(4.11) 



(4.12) 



Furthermore, the following Q-independent bounds will be established in Sec. |^. 
Lemma 4.2 For any Q G Had {J^, T), 

<i,^(A,A') <X%{X,X') V(A,y)GM^ 



(4.13) 



where Xab is independent ofQ and is defined in terms of the reference two-point function 
Up by 



xV^(A,A')=Fj(A)F^(A') 
where Y^{X) is the positive square root of 



72Maa^^p 



(—A, A) (no sum on A). 



(4.14) 



(4.15) 



Furthermore, ^^(A) (resp. Yj(A)j decays rapidly as X ^ +00 (resp. X —00) and is of 
polynomially hounded growth as X ^ —00 (resp. X +00). 
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Remarks, (i) The right-hand side of Eq. ( |4.15| ) is nonnegative because uaa<-^p is of 
positive type as a scalar bi-distribution, and this property is inherited under pull-back by 
72 (see Theorem 2.2 in |^). In fact, the calculation 

UAB^pif" ® F) = 4(/ V ® if'uBy) = a;o(^(/ V)^(/''^b)*) > (4.16) 

(summing on A and B) for /"^ G C^{M) (A = 1, . . . ,4) shows that uab^p is of positive 
type as a matrix-valued distribution. Positive type of uaa^^p follows in consequence. A 
similar argument (using ( [4.10|) and the property Q >0 for Q G Had(^, F)) shows that 
uab^'^', uab^^^'^ and their pull-backs by 72 share the matrix positive type property. 

(ii) The statements on the growth of are obtained from the Paley- Wiener-Schwartz 



theorem |2J] which entails that the Fourier transform of a compactly supported distri- 
bution is of at worst polynomial growth. Below, we will frequently use the fact that 
the product of a rapidly decaying function and one of polynomial growth is itself rapidly 
decaying. 



Because the bounds obtained in Lemma exhibit different behaviour in the four 
quadrants Ci, . . . ,6*4 of the (A, A')-plane it is convenient to decompose the averaged energy 
density ( |4.4|) as X = ^X^, where X^ is the contribution arising from quadrant Ck- We 
proceed to bound the X^ in turn. 

Starting with the second and fourth quadrants C2 = x and C4 = M''" x ]R~, 
Lemma ^]2| yields the Q-independent bound 

W^^i?^(A,A')|<5^xt(A,A') (4.17) 

in which each summand on the right-hand side is of at worst polynomial growth. This may 
be combined with the rapid decay of J^^ away from the diagonal to yield the following 
Q-independent bound on the contribution from these quadrants: 

\12+^a\< f d\d\'\j^^{\\')\Y^X%{\\') <oo. (4.18) 

JC2UC4 

We are left with the first and third quadrants. Since J"^^ exhibits polynomial growth 
along the diagonal, the previous argument will not allow us to bound all the terms arising 
from the decomposition ( 4.11 ). To see this, note that Lemma TJ2 applied to the P^QP^ 



term W\'j[ gives a bound growing polynomiallyQ in all directions in the first quad- 
rant. Similarly, the bound X'^^ for the PQ^P term is polynomially growing in the third 
quadrant. However, Lemma [4.2| suffices to bound the other contributions to Xi and X3 



because at least one factor in the relevant X\^ is rapidly decaying. Thus 

|Xl-7^l|< / rfArfA'|j^^(A,A')| [X;^b(A,A')+X/^(A,A')+X^^(A,A')] < 00 (4.19) 
and 

|X3-7^3|< / rfAciA'|j^^(A,A')| [X^r(A,A')+X/^(A,A')+X^^(A,A')] <oo, 

(4.20) 



^Note that it is the bound which is polynomially growing; we expect that is actually decaying. 
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where the remaining terms are 

7^l = /" d\dX'J^^{\,\')Wlj[{\,\') (4.21) 



and 

7^3 = -/ c/AdA'J^^(A,A')l^AB(A,A') (4.22) 

(the leading minus sign arises because it is —PQ^P which appears in (|4.8| )). The quantities 
TZi and 71^ are real, because J^'^, and W)[^ are hermitian matrix kernels. 

To complete the proof of the QWEI it is required to show that TZi and TZ^ are bounded 
from below independently of Q. We will present the argument for TZi in detail and indicate 
how the proof is modified for TZ^. Let xa (A G M"*") be a family of smooth, real- valued, 
nonincreasing functions such that xa equals unity on [0,A] and vanishes on [A + l,oo). 
It is clear that 

7^l=lim / dAdA'xA(A)xA(A')J^^(A,A')IVi/(A,A') (4.23) 
Jci 

and the cut-off xa now allows us to interpret the integral as a trace in the following way. 
Define a : M+ ^ M by 

a{\) = {l + \y/'[l + \Y^{\)y], (4.24) 

where | ■ |c4 is the usual vector norm on C'^. Then a is smooth, positive, bounded away 
from zero and of polynomially bounded growth. Next, define operators Ja and W on 
L2(M+,rfA) by 

^ /"OO 

(Ja^)^(A)= / dA'aA(A)(7A(A')J^^(A,A')W''(A') (4.25) 

and 







Jo cr{A')a{A") 

where have written o"a(A) = xa(A)o"(A). Now Ja is Hilbert-Schmidt (due to the cut-off) 
and self-adjoint while the properties of W are summarised in the following proposition, 
which is proved at the end of this section. 

Proposition 4.3 For all Q G Had(Jf , F), W is a positive trace-class operator with 

< Tr < ^ . (4.27) 



We may therefore rewrite Eq. (|4.23|) in the form 



7^l = lim Tr JaW . (4.28) 

A^oo 
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Introducing an orthornormal basis of eigenvectors Vn for W, and using W > and ( [4.27| ), 
we have 

Tr JaW = ^{Vn\ JkWVn) = ^{Vn \ J^Vn) K | Wv^) 
n n 

> inf spec ( Ja) ^ {vn \ Wvn) 

n 

> inf spec ( Ja) Tr W 

> ^ min{0, inf spec (Ja)}. (4.29) 

Noting that the right-hand side has no Q-dependence, the required lower bound on TZi 
now follows from the following proposition, which is proved in Sec. ^ 

Proposition 4.4 The spectrum of is bounded from below uniformly in A (with a finite 
lower bound). 

Turning to the integral TZ^, we follow exactly the same argument but with the single 
difference that the kernel J^^ on x M~ defines an operator on L^(]R~) <^C^ which may 
be bounded above by a nonnegative quantity; this is compensated by the leading minus 
sign in the definition of 71^. Thus TZi + TZs has a finite Q- independent lower bound as 



required and the proof of Theorem 4. 1 is complete. 



Proof of Proposition \^.,% - Using Lemma |4.2| and Cauchy-Schwarz, we first estimate 

'^^^^^''^'^^ - (1 + x'^l^'!t^x')y] I''" (1 + x''^ylt+ Ir^xny] 

(4.30) 

for ip G L^(R"'") ® C^, from which we obtain HW^v^H < |7r||(y9||. Thus W is bounded. A 
short calculation shows that 

{<f I W<f) = {i^ucBUj'^ ■ (V^^ ^) > (4.31) 

for G C^(M+) ® where 



1 

-if 
a 



r 



G yiR) ® , (4.32) 



and we have used the matrix positive type property and compact support of 72Ucb'^^ 
Thus W is positive. In order to show that W is trace-class it is now enough (by the 
lemma^ following Theorem XI. 31 in [^) to show that the formal trace of W is finite 
(whereupon the formal trace is indeed the trace). But this is just 

' dX a{Xr'6^^WBc{X, X) < r dX , .jr^^ill'^M < I ' (4-33) 



~'0 



'i + x^){i + \Y^{x)yY - 2 

where we have used Lemma |4.2| and thejfact that the diagonal of the kernel of a positive 

trace-class operator is positive. Thus W is trace-class on //^(M"*") ® C"^; putting this 

together with the positivity property, we have 0<Triy<7r/2as required. ■ 

*For this purpose, we regard L^(M+, dA) ® as L?{X, d\ ® d^j) where X is the locahy compact space 
M+ X Z4 and fx is the counting measure on Z4. The measure d\ ^ d/Lt is a Baire measure and the lemma 
may be appHed. 
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5 Proof of Lemma |^t2 



We first establish the existence of the quantities and ^^(A) defined by (^4.12 



and ( [4.151 ). Using self-adjointness of P", Cauchy-Schwarz and \\Q\\ < 1 (following from 
property (I) of Sec. |2]^ for Q G Had {J^, F)) we obtain from ([4.101) the inequality 



\J^\h^ft)\'< 



PT 




/i 







for any G 



1,2). This inequality underlies the following lemma, which will be 



proved at the end of this section. 

Lemma 5.1 The wave-front set of to'^^^ satisfies 



WF (J^^) C [j {p}x Ml X {p'} X -Ml, 

p,p'<^M 



(5.2) 



as a subset ofT*{M x M). 

This is by no means a sharp estimate of the wave-front set, but it will suffice for our 
purposes. 

Defining uab as in (|4.2|), uab^'^'^^ is a scalar bi-distribution with wave-front set con- 



tained in the right-hand-side of Eq. (|5.2|) . Now by Theorem 2.5.11' in the pull-back 
72My4_Ba;'' ^ exists provided the intersection of its wave-front set with the set of normals 
iV^2 of 72 is empty. One may show that 



72 



{(7(r),e;7(r'),n G T\M x M) | ea«"(r) = ^'yu'' {r') = 0} 



(5.3) 



(see Sec. 3 of |^ where a corresponding argument is given). This has trivial intersection 
with WF (uyisu;'*^''), because no null covector can annihilate a timelike vector. Accord- 
ingly, 72Ma_bi^^ exists in ^'(R^), and (again by Theorem 2.5.11' in |2^) its wave-front 
set obeys 



WF {'^luABOJ 



C 



X K X 



(5.4) 



Since uab is compactly supported, we may take Fourier transforms and conclude that 
the W]^Q do indeed exist. Exactly the same argument, using the microlocal spectrum 
condition ( p.50| ) in place of (p.2|), shows that the pull-backs ^2'^ab^\> exist with 



WF (72*mabc^p) C B 
It remains to prove Lemmas [4.2| and p.l| . 



X K X 



(5.5) 



Proof of Lemma K-^ ' An argument using regularising sequences in analogy with the proof 
of Theorem 2.2 in P] shows that the inequality (|5.1| ) is inherited by the pull-back and 
becomes 



\j;UABCO^^\fl ® /2)r < 72*«Aa4(/i ® /l) l*2UBBC0'p{f2 ® /2) V /i, G 



(5.6) 
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(no sum on either A or B). Substituting /i(t) = e~**'^, /2(t') = e~**''^', the required 
bounds (|4.13|) are obtained. 

As 'y2'^AA<^p is compactly supported, its set of singular directions (those directions in 
which its Fourier transform fails to decay rapidly) is given by 

S(7*MAA^p) = {(A, A') I (r, A; r', A') G WF {-f;uAAUJp) for some (r, r') G M^} 

= M+ X (5.7) 

(see Proposition 8.1.3 in |2^). Thus (—1, 1) is not a singular direction for 'J2Uaa<^p and we 
deduce that ^^(A) is rapidly decaying at A — +oo and of polynomially bounded growth 
as A — s> — oo by the Paley- Wiener-Schwartz theorem p4|. An analogous argument shows 
that yJ(A) decays rapidly as A ^ — oo and is polynomially bounded as A — +oo. ■ 
Proof of Lemma \5. 1\ : Suppose {pi,^i',P2, —^2) ^ WF (u;'"'''). We will show that 



{pi,^i;Pi,-^i)eWF{Jp) and (p2,6;P2,-6) e WF(4), (5.8) 

from which the required result now follows by the microlocal spectrum condition ( 2.501 ). To 
prove ( |5l8| ) fix charts (f/j, Kj) with pi G f/, {i = 1,2), and define ki so that = ^K[{pi)ki. 
Let Xi be arbitrary smooth spinor fields compactly supported in Ui with Xi{Pi) 0- 
By definition of the wave-front set on manifolds, (fci,— /C2) is a singular direction for 
Xi2i^'^^^ o K12, where we use the notation 

Xij = Xi® X^] f^ij = /ti X Kj . (5.9) 
Thus, if Vj are arbitrary conical neighbourhoods of kj, there must exist k'j G Vj such that 

[xuUJ^^'oK^^Ti^K^'C^k'^) (5.10) 
is not of rapid decay as a — * +00. Applying ( |0| ) to 

= f _, X.e-(-^"^'M ° , (5.11) 



and recalling (|1.3|) , we obtain 

[xi2UJ^^^ o f^u]^ iak[, -ak'2) < xu^poi^ii (aK, -'^K) [X22^p o ^22]^ {ak^, -ak^) 

(5.12) 

and since the left-hand side is not of rapid decrease as a — 00 we conclude that the same 
must be true of both factors on the right-hand side. Hence (as the Vj were arbitrary) 
(fci, —ki) is a singular direction for Xii^p ° ^^u and {k2, — ^^2) is a singular direction for 
X22^p ° 1^22 ■ Letting the support of each Xi shrink to {pi}, we have 

(/€(pi), k,; k{p,), -k,) G WF {Jp o «:-/) (5.13) 

and 

(«:(P2), k2; k{p2), -k2) G WF {Jp o /t^s^) , (5.14) 
from which (p.8[) follows immediately. ■ 
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6 Proof of Proposition |^74 



We now prove that the operators J\ on L^(M+, dX)®'C^ are bounded from below uniformly 
in A. To begin, we consider the related operator Ki^ acting on L^(]R"^, dX) by 



oo 



{K^^){\)= I rfA^^^^^A(^(A + A')/(A-A'MA'). (6.1) 







If the spin-connection terms vanished, J\ would be equal to ® 1. Our analysis of K\ 
is based on the following identity 

Lemma 6.1 If f = for real-valued g G C^(]R), then 

(A + A')/(A - A') = r ^/i?(A - /x)?(A' - /z) . (6.2) 

J-oo 

Proof: Note first that the right-hand side (RHS) exists for each A, A' G M. Changing 
variables to z/ = /i — (A + A')/2 and writing ( = (A — A')/2 

RHS of (6.2) = r ^ f + u] g{( - u)g{( + u) 



„ _ TT V 2 

(A + A')(?^?)(2C) 

(A + A')/(A - A') (6.3) 



as required, where we have also used g{u) = g{—u), the fact that i'g{C ~ ^)diC + ^) is odd, 
and a further change of variables. In addition, we have used -k to denote the convolution 
(/ii*/i2)(A) = / dX'/{2n)hi{X- X')h2{X'). m 

It follows from this identity that the kernel of K/^ may be rewritten in the form 

i^A(A,A') = ^^^^^^^3^^^^^ r dfx ixg{X - ix)g{X' - fi) . (6.4) 

We now define to have the kernels 

(A, A') = ^^^^.^^3^^^^^ r d^^ \^,\giX - /i)?(A' - /.) (6.5) 



and 



Kxix,x') = -^-^^^^^^ f d^.^^j{X-^.)W^) 



(6.6) 



The integrals in these kernels are bounded on compact subsets of M"*" x M"*", so the cut-off 
functions cta ensure that K^^ and are Hilbert-Schmidt. Clearly K\ = K\ — K~^] 
furthermore, the easily proven identity 



{if \K^(f) = dfi — dX'g{X' + fi)aA{X')^{X') 



(6.7) 
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oo / 

\ 



')) shows that Kj^ is positive. A similar argument estabhshes 



(vahd, say, for ip E C'^ 
positivity of K'l. 

Our aim is now to find a bound on K]^ which will allow a bound uniform in A to be 
obtained. (The operator K'^ is bounded for each A, e.g., by its Hilbert-Schmidt norm, 
but becomes unbounded in the limit A — > oo). Regarding the inner integral in ( |6.7| ) as an 
L^-inner product and applying Cauchy-Schwarz, we obtain 



(6. 



where 



47r3 



dfidX' fi ^(/i + A')cta(A') 



R+XR+ 



du\g{u)\^FA{u) 



and 



Fa(u) 



dX'{u - A')aA(A') 



/\2 



(6.9) 



(6.10) 



is bounded and nonnegative. Let us observe that this step depends in an essential way on 
the fact that, for > 0, the argument of 'g in (|6.7|) is bounded away from zero, together 
with the rapid decay property of 'g. 

The above analysis entails that —C\ is a lower bound for K\, but this is certainly not 
the sharpest bound. In fact essentially the same argument applies if Kj^ is replaced by 
and is adjusted to maintain K\ = — K^, with the conclusion that —C\/2 
is also a lower bound for Kj>^. The convenience of the choices made above is that, as we 
now show, the operator = Ja — K\ (g) U is form bounded relative to with relative 
bound no greater than i. To this end, we first use the convolution theorem to write 



1 r°° 

(y, I Lav) = 2 y dr[aAvY {r)^ f{r)e{r)[aAvY {r) 



for if G C^(]R^) eg) C^, where f denotes the matrix hermitian conjugate, 
sup II 6'(r) II c4, we then estimate 



(6.11) 
Setting C = 



\{v\Lav)\ < - / dT\{g[aAvY){T)\t,, 



< — 



2 
2 

C 

T 



^l(^7M")"(^)lc4 



d^ 
2^ 



dX_ 
2^ 



g{^i - A')(Ta(A')v'(A') 



(6.12) 



C4 



By comparison with the definition of K^, this implies 

d^ 



\{V I Lav)\ < \{V I {Kl®l)v)+ ^ 



dX'g{ix - AVa(A')v^(A') 



(6.13) 
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Using a similar argument to that used to obtain (|6.8|), the last term may be bounded by 
Cy^||(^|p where 



p poo 

du\g{u)\^Gf,{u) (6.14) 



c 



and the bounded nonnegative function Ga is given on [— C", oo) by 



C 



U+C 



Ga{u) = -— dX'aj,{X'y. (6.15) 



max{u— C,0} 



Thus La is form bounded relative to Kj[ as claimed above. Since J\ = — + La, 
we have 

\{ip\{J^~Kl®l)ip)\ < \{ip\L^ip)\ + \{ip\Klip)\ 

< ]^{^\{Kl®l)^) + {C^ + Cl)\M' (6.16) 

and, as is positive, it follows that Ja is bounded below by — (Ca + Ca). 

To complete the proof of Proposition [4.4| we must show that Ca and may both 
be bounded above uniformly in A. This holds because cta(A) is pointwise dominated by 
(j(A) and thus Fa and Ga are pointwise dominated by the functions F and G obtained 
by replacing cta by a in ( |6.10 ) and ( |6.15| ). Furthermore, F and G are of polynomially 
bounded growth so we obtain the bounds 



Ck< du\g{u)Y F {u) <oo (6.17) 
Jq 

and 

/oo 
du\g{u)\^ G(,u) < oo (6.18) 
■C" 



-c 

for all A, where the rapid decay property of g has been used. 



7 Majorana fields 

In this section we will indicate the changes required in Sec. ^ when treating Majorana 
fields. Suppose that 7o, . . . ,73 belong to a Majorana representation and that is the 
completion of ^cosp/ker 5'<] with scalar product given by (/i, h')s = {S^h, h'), cf. Remark 
(iii) at the end of Sec. pT2| . The field operators are then given by 

^(h) = B{q{h))- m+{f) = m{f+y = M/(r/+) h e ^^osp, / e %^ (7.i) 

where q : ^cosp — > -^cosp/ker S<i is the quotient map. Now, as in Sec. || one may choose, in 
a tubular neighbourhood of any timelike curve 7, induced frames (eo, ... ,63) and {Ea) 



25 



so that Co 1 7 = M, the tangent of 7 in proper time parametrization. However, it is now 
convenient to set Va = iE^, for then we have the two properties 

d'^^VA O = 7o and Tv^ = v\ (7.2) 

as we are working in a Majorana representation. 

To quantize the point-spht energy density (|3.4|) , we may substitute from ( |7.1|) and use 
the formula r[V ■ (eoWA/)]+ = V ■ (feov^) for / G ^(M) to obtain 

T(/ ®g) = 5^^^ (vI/(V ■ [feov\])^{vt9) - ■ [eov^)) (7.3) 

as the replacement for ( p.6|) . The CAR's may be used to show that (: T :)uj{f ® g) is 

symmetric in f,g. 

Using the formula Cq ■ Vt;^ = ao^f J, it follows that {T)^ may be expressed in the form 

{T)^ = C^''v+0v+uj2, (7.4) 
where uJ2{h, h') = u;(^(/i)^(/i')) is the two-point function and 

Z:^^ = i (1 ® «eo ■ V - zeo ■ V ® U) + -6^^ , (7.5) 
2 2 

with 

qAb _ i[S^^^ ® 1 - 1 ® 5^^^] . (7.6) 

The components of ctq are real in a Majorana representation and it follows that O^^ij) = 
720'^'^(t, r) is hermitian for each r. Moreover, ® "^^(^2 is of matrix positive type as 
shown by the calculation 

{vi ® 4^2)(F ® n = ^2(rK/^] ® t;+/^) = uj{<i^{vify^{v+f^)) > o (7.7) 

in which we have used Tv^ = v'^ and ( [7.1|) . 

From this point onwards, one may proceed to define (: p ■.)uj as in Sec. 3, and the 
statement and proof of Thm. 4.1 carry over (apart from some obvious changes) upon 
observing that Hadamard states u of the Majorana field obey 

WF(u;2) = {{p,^;p',-a e f*{M X M) I (p,0 ~ (p',0; e e AT+j (7.8) 

(cf. ||3^, note again that in this reference a different sign convention for the Fourier 
transform was chosen which results in the opposite form of the wave-front set there). 

8 Conclusion 

In this paper, we have established general QWEIs for the Dirac and Majorana fields in 
globally hyperbolic spacetimes. We conclude with various remarks. First, these QWEIs 
hold despite the fact that the 'classical' Dirac equation fails to obey the weak energy 
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condition. This is encouraging evidence that QWEIs are a widespread feature of all 
quantum field theories and that they are the correct replacement for the classical energy 
conditions. It would be interesting to understand whether QWEIs can be obtained in a 
general axiomatic setting. 

A second point was posed to us by Buchholz (private communication): given a weight 
in if the corresponding averaged energy density may be defined as a symmetric operator 
on a suitable dense domain (such as the domain of microlocal smoothness introduced in [§) 
in some Hilbert space representation. The force of our result (and the corresponding result 
in IP) is that such operators are semibounded, and therefore admit self-adjoint extensions 
(in particular, the Friedrichs extension). Can one give any interpretation to the evolution 
generated by this operator? The answer is not clear, but we speculate that there could 
be a connection with the dynamics discussed by Keyl in his recent study of quantum 



fields on timelike curves. The connection is somewhat tentative (in particular, the role 
of the weight must be understood), and would only be expected to hold under restricted 
conditions such as for static trajectories in static spacetimes. Nonetheless, it remains an 
intriguing possibility. 

Finally, although our approach does lead to explicit lower bounds on the various 
contributions to the averaged energy density, these are not expected to be optimal. We 
hope to return to this question elsewhere. 
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